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Abstract 

We consider = 1 supersymmetric Yang-Mills theory with fundamental matter in the large- 
Nc approximation in 1+1 dimensions. We add a Chern-Simons term to give the adjoint partons 
a mass and solve for the meson bound states. Here "mesons" are color-singlet states with two 
partons in the fundamental representation. The spectrum is exactly supersymmetric, and there is 
complete degeneracy between the fermion and boson bound states. We find that the mass spectrum 
is composed of two distinct bands. We analyze the properties of the bound states in each band 
and find a number of interesting properties of these states. In both bands, some of the states are 
nearly pure quark-gluon bound states while others are nearly pure squark-gluon bound states. The 
structure functions of many of the bound states found are very strongly peaked near x = 0. The 
convergence of the numerical approximation appears to be very good in all cases. 
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I. INTRODUCTION 



At this point there is no direct evidence for supersymmetry. Nevertheless supersymmetry 
is such a beautiful symmetry and provides such elegant solutions to a host of theoretical 
problems that many believe that it must be present in nature. It is a pressing experimental 
issue to see if nature takes advantage of this elegant option. Of course, we already know 
that supersymmetry is rather badly broken, since we do not see any superpartners for the 
particles of the standard model. It is assumed that all of the superpartners are in fact very 
heavy and that we will see them as we go to higher energies in accelerators. It is therefore 
extremely interesting to investigate what the properties of supersymmetric bound states 
might look like. There are indications that some unusual things can happen in a theory 
with supersymmetry 0, 1^ . 

We present here the solution of an A/" = 1 super Yang-Mills (SYM) theory with fun- 
damental matter in 1+1 dimensions. We also add a Chern-Simons (CS) term to give the 
adjoint matter a mass. There are several motivations to look at (l+l)-dimensional models, 
't Hooft [3| showed long ago that two-dimensional models can be powerful laboratories for 
the study of the bound-state problem. These models remain popular to this day because 
they are easy to solve and share many properties with their four-dimensional cousins, most 
notably stable bound states. Supersymmetric two-dimensional models are particularly at- 
tractive, since they are also super-renormalizable. Given that the dynamics of the gauge 
field is responsible for the strong interaction and for the formation of bound states, it comes 
as no surprise that a great deal of effort has gone into the investigation of bound states 
of pure glue in supersymmetric models 0]. Extensive study of the meson spectrum of 
non- supersymmetric theories has been done (see p] for a review). Recently an initial study 
addressed some of these states in the context of supersymmetric models [a]. In addition to 
these reasons for studying (l+l)-dimensional models, we are also interested in investigating 
them as the first step toward solving (2+l)-dimensional models. A (2+l)-dimensional model 
can be much more realistic, since it includes some transverse dynamics. 

Currently there are several numerical approaches to solving field theories. For QCD- 
like theories, lattice gauge theory is probably the most popular approach, since the lattice 
approximation does not break the most important symmetry, gauge symmetry. Similarly 
for supersymmetric theories, supersymmetric discrete light-cone quantization (SDLCQ) p, 
0, H, [l3| is probably the most powerful approach since the discretization does not break the 



most important symmetry, supersymmetry [ll|. In this paper we consider supersymmetric 



theories and follow this latter approach. To simplify the calculation we will consider only 
the laige-Nc limit j^, which has proven to be a powerful approximation for bound-state 
calculations. While baryons can be constructed in this limit they have an infinite 



number of partons, and thus practical calculations for such states are complicated. 

Throughout this paper we use the word "meson" to indicate the group structure of the 
state. Namely, we define a meson as a bound state whose wave function can be written as a 
linear combination of parton chains, each chain starting and ending with a creation operator 
in the fundamental representation. In supersymmetric theories, the states defined this way 
can have either bosonic or fermionic statistics. 

Previously we saw that the lightest bound states in A/" = 1 supersymmetric theories are 
very interesting In SYM theories in 1+1 and 2+1 dimensions, the lightest bound states 

in the spectrum are massless Bogomol'nyi-Prasad-Sommerfield (BPS) bound states j^. 
These states are exactly massless at all couplings. When we add a CS term to the (1+1)- 
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dimensional SYM theory, which gives a mass to the constituents, we find approximate BPS 
states. The masses of these states are approximately independent of the coupling, and 
at strong coupling these states are the lightest bound states in the theory 0. In (2+1)- 
dimensional SYM-CS theory we also find that at strong coupling there are anomalously light 
bound states [3]. In both 2+1 and 1+1 dimensions, these interesting states appear because 
of the exact BPS symmetry in the underlying SYM theory. We have also looked at the 
lightest bound states of SYM-CS theory with fundamental matter 3|- Again we see that 



the lightest bound states are significantly lighter than one would have naively expected. We 
found that the lightest state is nearly massless and well below threshold. 

We will see that this model has a number of interesting bound states in addition to the 
lowest mass state that we studied previously The bound states separate into two bands, 
a low-mass band and a high-mass band. Interestingly the low mass permits two solutions for 
bound states. The preferred solution has an unusual oscillatory convergence. Some of the 
states in the low-mass band are very light and well below threshold. The upper band has the 
standard linear convergence. In both bands, some of the states are nearly pure quark-gluon 
states, and some are nearly pure squark and gluons. Some have structure functions that are 
very sharply peaked at small longitudinal momentum fractions, and some have several of 
these properties. 

Throughout this paper we completely ignore the zero-mode problem jl3 . 15l |: however, it 



is clear that considerable progress on this issue could be made following our earlier work on 
the zero modes of the two-dimensional supersymmetric model with only adjoint fields [l6l |. 

The paper has the following organization. In Sec. Ill A) we consider three-dimensional su- 
persymmetric QCD (SQCD) with a CS term and dimensionally reduce it to 1+1 dimensions. 
We perform the light-cone quantization of the resulting theory by applying canonical com- 
mutation relations at fixed x"*" = (a;*^ + x^)/a/2 and choosing the light-cone gauge {A~^ = 0) 
for the vector field. After solving the constraint equations, we obtain a model containing four 
dynamical fields. We construct the supercharge for this dimensionally reduced theory. In 
Sec. Ill BI we discuss the structure of the lighter meson bound states in the large-iV^ approx- 



imation. In Sec. Ill CI we discuss the addition of a CS term to the supercharge IT], |l8|, uS 
We explain that, in the context of this model, this is equivalent to adding a mass to the 
partons in the adjoint representation. In Sec. II I II we discuss the bound-state solutions that 
we find, including the mass spectrum, structure functions and the numerical convergence of 
the various bound states. In Sec. |TV| we discuss our results and the future directions that 
are indicated by this research. 



II. SUPERSYMMETRIC SYSTEMS WITH FUNDAMENTAL MATTER 
A. Construction of the supercharge 

We consider the supersymmetric models in two dimensions which can be obtained as the 
result of dimensional reduction of SQCD2+1. Our starting point is the three-dimensional 
action 

S = J d^xtT (^~F^^F>"' + ^AVD^A + D^^^D^'i + i^D^V^'■<^ 

-g[^Ki + eA-^]) . (2.1) 
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This action describes a system of a gauge field A^, representing gluons, and its superpartner 
A, representing gluinos, both taking values in the adjoint representation of SU{Nc), and two 
complex fields, a scalar ^ representing squarks and a Dirac fermion \1/ representing quarks, 
transforming according to the fundamental representation of the same group. In matrix 
notation the covariant derivatives are given by 

D^A = d^A + tg[A^, A] , D^^ = d^^ + tgA^^ , D^^ = d^^ + tgA^^ . (2.2) 

The action ()2.H) is invariant under the following supersymmetry transformations, which are 
parameterized by a two-component Majorana fermion e: 

SA^='-eT,A, 5A = ^F^,r'^^£, 

S^='-e^, <5vi> = -lr%D^e. (2.3) 

Using standard techniques one can construct the Noether current corresponding to these 
transformations as 

+ ^D^^T^"'e^ . (2.4) 

We will consider the reduction of this system to two dimensions, which means that the 
field configurations are assumed to be independent of the transverse coordinate x^. In 
the resulting two-dimensional system we will implement light-cone quantization, where the 
initial conditions as well as canonical commutation relations are imposed on a light-like 
surface = const. In particular, we construct the supercharge by integrating the current 
()2.4|) over the light-like surface to obtain 

- '-^eD+^ + ^^D.^T+'^e^^ . (2.5) 

Since all fields are assumed to be independent of x^, the integration over this coordinate 
gives just a constant factor, which we absorb by a field redefinition. 

If we use the following specific representation for the Dirac matrices in three dimensions: 

r° = a2, r' = tai, r^ = tcT3, (2.6) 

the Majorana fermion A can be chosen to be real. It is also convenient to write the fermion 
fields and the supercharge in component form as 

A=(a,a)^, ^=(i:,tpy, Q={Q^,Q-f. (2.7) 
In terms of this decomposition the superalgebra has an explicit (1, 1) form 

{Q+,Q+} = 2V2P+ , {Q-,Q-} = 2V2P- , {Q+,Q-} = 0. (2.8) 
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The SDLCQ method exploits this superalgebra by constructing P~ from a discrete approx- 
imation to Q~ rather than directly discretizing P~, as is done in ordinary DLCQ 0]. 

To begin to eliminate nondynamical fields, we impose the light-cone gauge (A"*" = 0). 
Then the supercharges are given by 

g+ = 2 1 dx- I^Xd^A' + - - '^^d^ij + '-d.^P^^^ , (2.9) 

Q- = -2 j dx- (^-A9_A- + ze^D2^-^^2^^e + ^5-(v^^e-e^^)) • 

(2.10) 

Note that apart from a total derivative these expressions involve only left-moving compo- 
nents of the fermions (A and ip). In fact, in the light-cone formulation only these components 
are dynamical. To see this we consider the equations of motion that follow from the action 
fl2.1|) . in light-cone gauge. Three of them serve as constraints rather than as dynamical 
statements; they are 

dSX = A] + ^e^t _ ,^^t) , (2.11) 



and 
with 



dlA-=gJ, (2.13) 

J = i[A\ d^A"^] + ^{A, A} - ihd-^^^ + i^d^^^ + v^^j^j^ . (2.14) 
v2 

Apart from the zero-mode problem one can invert the last constraint to write the 
auxiliary field A~ in terms of physical degrees of freedom. 

In order to solve the bound-state problem 2P~^P~\M) = M^|M), we apply the methods 
of SDLCQ. Namely we compactify the two-dimensional theory on a light-like circle (— L < 
x^ < L), and impose periodic boundary conditions on all physical fields. This leads to the 
following mode expansions: 

1 oo ^ 

4(0, X-) = E 7^ (a.,(A:)e-^'='^-/^ + a],ik)e^'-^- , (2.15) 

k=l 
^ oo 

A.,(0,x-) = -^-==$^(%(A;)e-'=--/^ + 6],(A;)e^'=--/^) , (2.16) 

m X-) = -= 5^ -= (Q(fc)e-^--"/^ + gl(A:)e^^-^-/^) , (2.17) 

k=l 
^ oo 

MO,x-) = , J2 {d^{k)e-'^^^''^ + d\{k)e'^^''''^^ . (2.18) 
2 4 V 2Z/ j^—Y 



We drop the zero modes of the fields; including them could lead to new and interesting 
effects (see (lij, for example), but this is beyond the scope of this work. 
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In the light-cone formahsm one treats as the time direction, thus the commutation 
relations between fields and their momenta are imposed on the surface = 0. For the 
system under consideration this means that 



[Ali0,x-),d.AU0,y-)] = I (^5u5kj - ^5,,5k^ 5{x- - y ) , 

[^i(0,x"),9_^j(0,?/")] = i5ij5{x' - y') , 
{V^,(0,x-),V^,(0,i/-)} = V25i,5{x- -y-). 

These relations can be rewritten in terms of creation and annihilation operators as 



6. 



13 ' 



Ci, C^j 



5. 



Si. 



dill • dj ■ 



6ij. 



(2.19) 

(2.20) 

(2.21) 
(2.22) 



(2.23) 



(2.24) 



In this paper we will discuss numerical results obtained in the large-A^^c limit, i.e. we ne- 
glect 1/Nc terms in the above expressions. Although 1/Nc corrections may have interesting 
consequences, they are beyond the scope of this work. 

Substituting the result into the expression for the supercharge and omitting the boundary 
term, we get 

Q- = Q- + gr + Q2 + • (2-25) 



where is the supercharge of pure adjoint matter [9[. The three supercharges that govern 
the behavior of the fundamental matter in these states are 



Q- = -A^j dx- (iV2id^e - ^V2^^^e) ^\ , 
Q^- = -^Jdx-{2^^^)^X, 



(2.26) 
(2.27) 
(2.28) 



After substituting the expansions ()2.15ll2.lHjl . one gets the mode decomposition of the su- 
percharge, 

z2-i/4^yi ^ 



TT 



ki,k2,ks=l 



(2.29) 
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with 



Q2 = 



+b],{k^)4{k2)c^{ks) - cl{k2)kj{k^)c^{ks)]Sk,M+k, , (2.30) 



dl{k2)blj{ki)dj{ks) + d]{ks)di{k2)bij{ki)]Sk, (2.31) 

^3" = Tr77=1=[^4(^i)ci(A;3)aij(A;2) +4(^3)4(A:2)rfj(A;i) + 



{d\{ki)aij{k2)cj{k'i) + aJj(A;2)c](A;3)(ii(A;i))5fci,it3+ik2 



{k'i)di{ki)aij{k2) + (i|(fci)a^(A;2)cj(A;3) + 
(cJ(/c3)ay(A;2)(ij(A;i) + at(/c2)(^](A;i)Q(A;3))(^fe3,iki+ik2] ■ (2-32) 

We have dropped from these expression the terms that connect the states that are closed 
loops of adjoint partons to the string-like meson states. These interactions are of order 
1 / y/^c and can be neglected in the large- A^c approximation. 



B. Meson structure 

We will consider here only meson-like states. In the large-iY^ approximation these are 
color-singlet states with exactly two partons in the fundamental representation. The boson 
bound states will have cither two bosons or two fermions in the fundamental representation. 
In general a boson bound state will have a combination of these types of contributions. Be- 
cause this theory and the numerical formalism are exactly supersymmetric, for each boson 
bound state there will be a degenerate bound state that is a fermion. The fermion bound 
state will have one fermion in the fundamental representation and one boson in the fun- 
damental representation. In the string interpretation of these theories, such states would 
correspond to open strings with freely moving cndpoints. In the language of QCD, the model 
corresponds to a system of interacting gluons and gluinos which bind dynamical (s) quarks 
and anti-(s)quarks. In the large- A^c limit we will have to consider only a single (s)quark- 
anti-(s)quark pair. Thus the Fock space is constructed from states of the following type: 

fl{ki)aii,{k2) . . .6L..(^n-i) . ..fliknm • (2.33) 

Here J} and f} each create one of the fundamental partons, and |0) is the vacuum annihilated 
by flij, bij, Ci, Ci, di, and di. 

The other color-singlet bound states in this theory are states that are composed of traces 
of only adjoint mesons. These can be considered to be loops. At finite Nc this theory has 
interactions that break these loops and insert a pair of fundamental partons, making an 
open-string state. This type of interaction can, of course, also form loops from open strings 
and break one open string into two open strings. In principle, a calculation of the spectrum 
of such a finite- A^c theory is within the reach of SDLCQ. The only significant change is to 
include states in the basis with more than one color trace. 
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C. Supersymmetric Chern Simons theory 



The CS term we use in this calculation is obtained by starting with a CS term in 2+1 
dimensions and reducing it to 1+1 dimensions. This term has the effect of adding a mass 
for the adjoint partons. In this calculation we are including fundamental matter because 
we are interested in QCD-like meson bound states. Without a mass for the adjoint matter, 
this theory is known to produce very long light chains of adjoint partons. In a finite-iVc 
calculation we would not have these very long chains because they would break, but in the 
large- A^'c approximation they do not. While SDLCQ can be used to do Unite-Nc calculations, 
it is much easier to add a mass to restrict the number of adjoint partons in our bound states. 
We choose the CS mechanism to give the adjoint partons a mass because we can do this 
without breaking the supersymmetry. 

The Lagrangian of this theory is 

^ = CsQCD + '^Ccs , (2.34) 
where jCsqcd is the SQCD Lagrangian we discussed earlier, k, is the CS coupling, and 

Ccs = e'^"' (^A^d^Ax + jgA^A.A^^ + 2^^^ . (2.35) 

A trace of the color matrices is understood. The constraint equation (j2.12j) gains a third 
term of the form — kA/a/2 on the right-hand side, and the definition of the current in (j2.14p 
now has an additional term, Kd-A"^. The discrete version of the CS part of the supercharge 
in 1+1 dimensions is 




(2.36) 



It is informative to compare Q^s with the one term in the supercharge for Af = 1 SYM 
in 2+1 dimensions [i^ which has an explicit dependence on the transverse momentum k±. 
This fc^-dependent term has the form 




(2.37) 



where k± = 2TTn±/L± is the discrete transverse momentum. Notice that k±_ and k enter the 
supercharge in very similar ways. Addition of a CS term has the effect of changing k± to 
k± + iK. Because the light-cone energy is of the form {kj_ + m'^)/k~^, k± behaves like a mass, 
and therefore n also behaves in many ways like a mass for the adjoint particles. 

The partons in the fundamental representation in this theory will remain massless. Of 
course, in a more physical theory the supersymmetry would be badly broken; the squark 
would acquire a large mass, and only the quarks would remain nearly massless. 

III. NUMERICAL RESULTS 

This SYM-CS theory with fundamental matter has two dimensionful parameters with di- 
mension of a mass squared, the YM coupling squared g^N^/ n and the CS coupling squared 
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FIG. 1: The mass-squared spectrum, with Z2 odd, in units of g'^Nc/'K, at a resolution of = 6 as 
a function of (a) k at g^/Wjir = 1 and (b) gy' Nc/tt at k = 1. 



k"^. The latter is also the mass squared of the partons in the adjoint representation. All 
of the masses in this paper will be given in units of gy^Nc/Tc, which will usually be sup- 
pressed. Furthermore, we are only considering meson bound states. These are states of the 
form shown in Eq. ()2.HHj) with two fundamental partons linked by partons in the adjoint 
representation. Since we are working in the large- A'c approximation, this class of states is 
disconnected from the other allowed class of pure adjoint matter bound states and multi- 
particle states. This theory also has a Z2 symmetry jij] which is very useful in labeling 
the states and reducing the dimension of the Fock basis that one has to consider in any one 
diagonalization step. For this theory the Z2 symmetry divides the basis into states with an 
even or odd number of gluons. 



A. Spectra 

We find that the spectra of meson bound states, for both odd and even Z2 symmetry, 
divides into two bands of states as we increase k, a light-mass band and a heavy-mass band, 
as can be seen in Fig. Ufa). As n grows, a gap develops in the spectrum. This mass gap 
reflects in part the number of massive adjoint partons in the bound states. At very large 
K, the low-mass band will comprise states with only fundamental partons. We will look at 
the bound states in both bands and consider various values of k, relative to gy^Nc/ir. In 
Fig. n^b) we see that at k = 1 this gap grows with the coupling. The lowest mass state 
remains very light even for values of g\/Nc/TT up to 2. As far as we can tell, this state 
remains very light even at very large couplings. 

In Fig. |21^a) and (b) we show the spectrum in the Z2-even sector, with g^fWjii = 1, as 
a function of resolution at k =0 and 1, respectively. We see that at large resolution with 
K = the two bands merge. It is unclear from this figure whether the states in the lower 
band mix with the states in the upper band in the region where their masses overlap. An 
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1/K 

(a) (b) 

FIG. 2: The mass-squared spectrum, with Z2 even, in units of g^Ncji^-, as a function of the 
resolution at (a) 5 = and (b) Nc/tt = 1 



inspection of adjacent states seems to show that these states remain on smooth trajectories 
and do not repel, as one would expect if they were interacting. This seems to indicate that 
separate analysis of these two sets of states is possible. For the most part, however, we will 
look at the situation where shown in Fig. |21^b). In this case the two bands are well 

separated even at the highest accessible resolution, and, therefore, there is not a problem 
considering the two bands separately. 

Let us first consider the bound states in the lower band. To understand the physics 
underlying these bound states we will look at both the spectrum and the structure functions 
of individual states. In this theory we have four species of partons: adjoint fermions and 
bosons, and fundamental fermions and bosons. Therefore, we will look at four structure 
functions that give the probability of finding a particular variety of parton at a particular 
value of longitudinal momentum fraction in a particular bound state. We need all of this 
information about each bound state to properly analyze this theory. 

Although the procedure that we follow for identifying bound states is the standard one 
in DLCQ, it is worthwhile to give a brief description of that procedure here because of 
the unusual nature of the lower band. To identify a bound state, we start at the lowest 
value of the resolution K, usually K = ?>, and select a particular eigenvalue, usually the 
lowest nonzero one. We then look at the properties of this state, which we always calculate 
along with the mass. In principle, we can look at the entire wave function, but for higher 
resolution there is more data than we can efficiently handle. It is efficient to look at the 
average number of partons of the various types and sometimes the average momentum of 
some of these partons. We then move to the next higher resolution and try to identify an 
eigenstate that has the same or very similar properties to the state at lower resolution. We 
continue this through all the resolutions. This then gives us the mass squared of a particular 
bound state and the various properties of the state as functions of the resolution. If this 
history of a state as a function of the resolution makes sense, and can be extrapolated to 
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K 




< n > 


< > 


< nfB > 


< naF > 


< Ufp > 


3 


0.178 


2.30 


0.30 


1.03 


0.01 


0.97 


4 


0.006 


2.56 


0.51 


1.86 


0.05 


0.14 


5 


0.049 


2.69 


0.63 


1.29 


0.06 


0.71 


6 


0.016 


2.83 


0.75 


1.71 


0.08 


0.30 


7 


0.029 


2.84 


0.76 


1.45 


0.08 


0.55 


8 


0.022 


2.92 


0.83 


1.58 


0.09 


0.42 


9 


0.025 


2.92 


0.83 


1.49 


0.10 


0.51 


10 


0.024 


2.96 


0.86 


1.52 


0.10 


0.48 


11 


0.025 


2.97 


0.87 


1.49 


0.11 


0.51 


12 


0.025 


3.00 


0.89 


1.48 


0.11 


0.52 


13 


0.026 


3.01 


0.90 


1.47 


0.11 


0.53 



TABLE I: Properties of the lowest mass boson bound state in the Z2-even sector, including the 
average numbers of adjoint bosons aB, adjoint fermions aF, fundamental bosons fB, and funda- 
mental fermions fF, for different values of the resolution K. The mass squared is given in 
units oi g^Nc/iT. The CS coupling is k = g^/Njir. 

infinite resolution, we say we have found a bound state. We then repeat the process to obtain 
another bound state. For the oscillatory states presented below, we actually started with the 
lightest state at the highest resolution and then followed the states down to low resolution. 
We then repeated the process with the next highest states at the highest resolution. 

The lower band of this model is unique relative to all other models we have studied 
in SDLCQ or DLCQ. We find that there are two different ways of carrying out the above 
procedure to find the states in the lower band. We can follow the eigenvalues from resolution 
to resolution in two ways. As a result we end up with two sets of bound states, all of whose 
properties are very smoothly behaved as functions of the resolution. Two such states and 
their properties are shown in Tables HI and ITTl The masses in Table H] are very small. We have 
discussed the lowest mass state in this spectrum in detail elsewhere 0]. The convergence 
as a function of the resolution is oscillatory, but clearly as we move to large resolution the 
convergence is very good. In Fig. El^a) we present the spectrum as a function of 1/ and 
we see that the oscillations lead to good convergence that extrapolates nicely to = cxd. 
These states have masses near zero. The lowest state has on average one adjoint parton; 
therefore, the threshold is at k^. Thus this is a deeply bound state. In Fig. Efb) we show the 
spectrum for the next highest mass bound state, and again we see very good convergence 
to a very hght state. In fact all of the data in the lower band describe a series of oscillatory 
states at higher masses. 

Shown in Table |n] is another way of organizing the same set of data used to analyze 
the oscillatory states discussed above, and all these states also converge numerically, but to 
infinite mass. As a function of the resolution we see almost no variation in the properties 
of this state. We will refer to these states as divergent states, and we will argue that they 
are divergent. This data and a fit to the data points are shown in Fig. |2^c), and there are 
no oscillations. In Fig. Efd) we see another state of this type. It appears that the divergent 
states do extrapolate to an infinite mass as — > oo. The curves, however, can be fit by a 
variety of functions. In fact the functions a + hj^/K ^ 1/ K and a + 61ogi^ +1/ K produce 
equally good fits. The fits shown are of the former form and have intercepts between 10 and 
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K 




<n> 


< naB > 


< njB > 


< naF > 


< Ufp > 


3 


0.29 


2.39 


0.33 


1.01 


0.06 


0.99 


4 


0.75 


2.58 


0.47 


0.91 


0.11 


1.09 


5 


1.23 


2.73 


0.59 


0.87 


0.14 


1.14 


6 


1.67 


2.82 


0.66 


0.85 


0.17 


1.15 


7 


2.06 


2.88 


0.70 


0.84 


0.18 


1.16 




2.41 


2.92 


0.72 


0.84 


0.20 


1.16 


9 


2.73 


2.94 


0.73 


0.84 


0.21 


1.16 


10 


3.02 


2.96 


0.74 


0.84 


0.22 


1.16 


11 


3.28 


2.97 


0.74 


0.84 


0.23 


1.16 


12 


3.53 


2.98 


0.74 


0.84 


0.23 


1.16 



TABLE II: Same as Tabled but for the lowest mass divergent state. 



11. In Fig. 13] we put all four plots from Fig. El on the same graph, so that one can clearly see 
that at low resolution both curves involve the same data. In the DLCQ procedure a state 
must extrapolate nicely to i^' = cxd, if it is a true bound state of the theory. It is unclear 
from the data whether this is true for these states. Our prejudice is that this bound state 
diverges as i^' — > cxd. 

In Fig. |S1 we show similar results for the Z2-odd sector. In Fig. IHl we again show the 
oscillatory and divergent fits on the same figure, but now we add a few more states to better 
show the interlacing of these fits to the same data. There are, of course, degenerate fermion 
bound states for all these states, whose properties we do not show. 

The oscillatory behavior made this calculation particularly challenging numerically. We 
were forced to go to very high resolution, K = 13, to be certain that the spectrum really 
converged. This was made even more difficult because we had four species of particles in 
the problem. 

To get some insight into the behavior of the low-mass band, we recall the behavior of 
the two-particle continuum in DLCQ. The "mass squared" of two free partons of mass m at 
resolution K is 

1 1 



K — n n 



(3.1) 



where one free parton has longitudinal momentum K — n and the other has n [22]. In DLCQ 
this formula produces a band of states. It is interesting to compare the top and bottom of 
this band with the top and bottom of the low-mass band of bound states we have been 
studying. The top of the band is obtained by fixing n at one. If we then take K large, we 
find = rn?{K +1 + 1/ K + ...). Not surprisingly, the divergent curves in Figs.Efc) and (d), 
and Figs.Efc) and (d) can be fit with a function of the form = a + bK + c/K. In essence, 
the form of the divergent state is comparable to the shape of the top of the continuum band. 
Similarly the bottom of the continuum band is obtained by taking n = K/2 for K even and 
n = {K — l)/2 for K odd. The bottom of the continuum band oscillates as a function of 
K. For K even it is 4m^, and for K odd it is 4m^(l + 1/K + ...). It again appears that 
the lowest mass states resemble the bottom of the continuum band. We would expect to 
see a connection between the weak-coupling bound states and the free spectrum, but this 
connection appears to extend to g'^Nc/n = k^. Based on this argument, the divergent states 
actually do diverge and are not true states of the theory. 
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(c) (d) 

FIG. 3: The mass squared in units of q^Nc/tt, as a function of 1/i^ for k = gyN^Ji: in the Z2-even 
sector, of (a) the lowest mass oscihatory state, (b) the second lowest mass oscillatory state, (c) the 
first divergent state, and (d) the second divergent state. The solid curve is a fit to the computed 
points. 

The lowest mass Z2-even state in the low-mass band has on average one adjoint parton, 
which has mass 1.0, and the fundamental partons are massless. One would therefore expect 
the lowest bound state in the spectrum to have a mass of order 1.0. In A/" = 1 SYM-CS 
theory 0, we previously found anomalously light states at strong coupling, but all of 
these states were at or close to threshold. Here the lowest mass state is anomalously light, 
in fact nearly massless, and therefore well below threshold. 

At yet stronger coupling, g'^Nc/n = IOk^, we find that oscillations are significantly 
stronger, as seen in Fig. [71 and at K = 12 the curve is not totally converged. The mass 
squared at K = 12 is 0.068, but without complete convergence we are not able to say any- 
thing precise about the variation of the bound-state mass as a function of the coupling from 
g^Njn = ^2 to g^Njn = IOk^. 

B. Structure functions 

The two largest structure functions for the lowest mass state are shown in Fig. |H| Even 
though this is an exactly supersymmetric theory, nearly all the adjoint partons are gluons. 
About 2/3 of the wave function of this state is composed of two fundamental bosons and an 
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K 

FIG. 4: The mass squared in units of g^Nc/n-, as a function of 1/K for k = q^JNc/tt in the Z2-even 
sector of oscihatory and divergent mass fits. 
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(c) (d) 
FIG. 5: Same as Fig. |3 but for the Z2-odd sector. 

adjoint parton, and about 1/3 of the wave function is made of two fundamental fermions. 
Within the context of the standard model, this state is primarily a bound state of two 
squarks and a gluon. We see that both of these distributions are peaked at small x. This 
reflects strong binding of the fundamental partons, allowing them to be widely separated 
in momentum, combined with only a small contribution to the momentum from the adjoint 
boson. 
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FIG. 6: Same as Fig. ^ but for the Z2-odd sector. 
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FIG. 7: The mass squared of the lowest mass oscihatory state in units of g^Nc/i:, as a function of 
1/K for IOk^ = g^Nc/TT in the Z2-even sector. 

The average number of partons in the next lowest state is about 4.5. For this state, 90% 
is composed of gluons and squarks, and the structure functions are shown in Fig. IHl The 
squarks have a significantly wider distribution than in the lowest state. The gluons remain 
very sharply peaked at small x. 

For the lowest state in the low-mass band with Z2 odd, the average number of partons 
is slightly less that 4. Again the partons in the adjoint representation are almost entirely 
bosons. However, the partons in the fundamental representation are almost entirely quarks, 
so this is the lightest standard-model meson. This is a deeply bound state, since the gluons 
have mass 1 and on average there are close to two gluons per bound state. The structure 
functions for the gluon and quarks are shown in Fig. ^1 As we have seen in all the states 
that we have considered, the gluons are strongly peaked at x near zero. The quarks have a 
strong peak at small x but have a significant distribution at other values of x. 

C. Spectra and structure functions in the upper band 

While the lower mass band appears to have a very unusual behavior as a function of the 
resolution, the bound states in the upper band behave very similarly to the bound states 
found in most SDLCQ calculations. The convergence is excellent, starting from the lowest 
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(a) (b) 

FIG. 8: Structure functions of the bound state = 0.025 at resolution K = 12, with Z2 even, for 

(a) adjoint bosons aB and (b) fundamental bosons with the CS coupling fixed at k = g\/ Nc/n. 
The solid curve is a fit to the computed points. 




1 . 2 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



(a) (b) 

FIG. 10: Structure functions of the bound state = 0.135 at resolution K = 12, with Z2 odd, 
for (a) fundamental bosons fB and (b) fundamental fermions fF, with the CS coupling fixed at 
K = g\J Nc/tt. The solid curve is a fit to the computed points. 
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(a) (b) 

FIG. 11: The mass squared in units of g'^Nc/ir, as a function of 1/K for k = g\J Nc/tt in the 
Z2-even sector, of (a) the second lowest mass state of the upper band, with = 9.03, and (b) 
the lowest mass state of the upper band, with = 8.31. The sohd curve is a fit to the computed 
points. 

resolution, and the data is fit very well as a function of l/f^ by a line with a small slope. 
In Fig. ^2 we show that the lowest two bound states in the upper band converge to masses 

= 8.31 and = 9.03 at infinite resolution. We carry the calculation out to resolution 
K = 13, but the structure functions are calculated at = 12. From the structure functions 
of the state at = 8.31, we find that this state is 85% quarks and gluons, so it is very 
much like a standard-model meson. The average number of partons is 3.98, so for the most 
part the bound state contains two gluons. It is very interesting that this bound state is 
very much like a QCD meson and is the solution of a supersymmetric field theory. The 
structure functions are shown in Fig. ^| While the gluons are peaked at small x, they have 
a much larger spread than we saw in the lower band. The fundamental fermions are very 
sharply peaked at small x. The structure functions of the state at = 9.03 are interesting 
because the state has nearly equal mixtures of all four constituents. The structure functions 
are shown in Fig. ^| The shapes are different, particularly for the distribution of adjoint 
fermions, which are spread over the region of larger x. 

In Fig. E]we see that in the Z2-odd sector we again have an excellent linear fit. The two 
lowest states have masses of = 8.93 and = 9.60. The structure functions for this 
sector are similar to those for the Z2-even states. The structure functions for the state with 

= 8.93 are shown in Fig. Upland are very similar to the structure functions for the state 
with = 9.03 in the Z2 -even sector. 

IV. DISCUSSION. 

In this paper we studied A/" = 1 SYM-CS theory with fundamental matter in 1+1 dimen- 
sions. The CS term was included to give masses to the adjoint partons. The calculations 
were performed at large A*",, in the framework of SDLCQ; namely, we compactified the light- 
like coordinate x~ on a finite circle and calculated the Hamiltonian as the square of a 
supercharge Q~, which we then diagonalized numerically. We found that the spectrum of 
this theory has two bands, a lower mass band and an upper mass band. With a CS term 
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(a) (b) 

FIG. 12: Structure functions of the bound state = 8.31, with Z2 even, for (a) adjoint bosons 
aB and (b) fundamental fermions /F, with the CS couphng fixed at k = g^J A^'c/vr. The sohd curve 
is a fit to the computed points. 

present, we found that these bands separate for g'^Nc/ir = k,^, and we can easily study them 
separately. 

For very massive adjoint partons, the lower mass band becomes a set of massless bound 
states composed of only fundamental partons. When the mass for the adjoint partons is 
reduced, the massless states become the states of the lower mass band. The states in the 
low-mass band are unusual in that their convergence is oscillatory. We show that the upper 
and lower bounds of the low-mass band have the same shape as the DLCQ approximation 
to the two-particle free spectrum, which has an oscillatory behavior for its lower bound and 
a growing upper bound. We argue therefore that the oscillatory behavior is a numerical 
remnant of the free two-particle spectrum. In previous work [13] we found that at g^Nc/ir = 

the lowest mass state of this theory was anomalously low, and in fact close to zero, while 
threshold for this state is at = 1. On average this state has one massive gluon and 
two squarks. This is the fifth supersymmetric theory where we have found anomalously 
light states. In SYM in 1-1-1 0] and 2+1 [2^ dimensions, there are massless BPS states. 
In SYM-CS theory in 1-1-1 dimensions we saw that at strong coupling the lightest states 
are approximately BPS states whose masses are independent of the YM coupling |2j. In 
SYM-CS theories in 2-1-1 dimensions at strong coupling, there is again an anomalously 
light bound state The next lightest state in this band has about four partons. We 
carry the calculation to resolution K = 13, so in principle these states could contain 13 
partons. In addition, almost all of the partons in the adjoint representation are gluons 
rather than gluinos. Also, the content of the particles in the fundamental representation 
is rather pure, either squarks or quarks. Thus some of the boson bound states are rather 
pure combinations of quarks and glue, in spite of the fact that they are bound states of an 
exactly supersymmetric theory. 

The bound states in the upper mass band are rather conventional SDLCQ bound states. 
They converge very quickly and have the conventional linear behavior in 1/K. Similar to the 
lower mass band, most of the partons in the adjoint representation are gluons. In addition, 
the states are again rather pure quark or squark states. When we take the adjoint parton 
mass large, these are the states that have at least one adjoint parton. 

We calculate the structure functions of the states in both bands for all four species of 
partons, and we find that most of the distributions are strongly peaked at x near zero. This 
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(c) (d) 

FIG. 13: Structure functions of the bound state = 9.03 for (a) adjoint bosons aB, (b) 
fundamental bosons fB, (c) adjoint fermions aF, and (d) fundamental fermions fF, with the CS 
couphng fixed at k = g\J Nc/i^ and Z2 even. The solid curve is a fit to the computed points. 

is an indication that for q'^Nc/t^ = these are strongly coupled bound states. In a weakly 
coupled bound state one would, on average, expect the partons to share the momentum 
fraction equaUy, and therefore the structure function would be peaked at 1/3 for a state 
with three partons, for example. Since several of the states we consider are almost pure 
quark-gluon bound states, they might not be affected by supersymmetry breaking, which 
would give large mass to the squarks and gluinos. It would be interesting to investigate in 
more detail the properties of these states and compare them with QCD. 

There remains a considerable amount of work to be done on SYM-CS theories with 
fundamental matter. The most straightforward extension of the present work is to consider 
calculations in 2+1 dimensions [20, 23, 24j. The M = 1 theory in 2+1 dimensions is easily 
within our reach. Beyond that the J\f = 2 theory ji^ in 2+1 dimensions, which is the 
dimensional reduction of the A/" = 1 theory in 3+1 dimensions, will be very interesting. For 
finite Nc we could also look at the mixing of closed-loop states and finite string-like meson 
states. For k = the bound states of all these theories provide an interesting field-theoretic 
model for strings. 
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FIG. 14: The mass squared in units of g'^Nc/ir, as a function oil/K for k = g^J Nc/tt in the Z2-odd 
sector, for (a) the lowest mass state of the upper band, with = 8.93, and (b) the second lowest 
mass state of the upper band, with = 9.6. The solid curve is a fit to the computed points. 




(c) (d) 
FIG. 15: Same as Fig. [13 but for = 8.93 and Z2 odd. 
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